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We reexamine the conventional arguments concerned with the scaling and decoupling solutions 
for the ghost and gluon propagators in the Landau gauge Yang-Mills theory. We point out a few 
issues to be clarified, which seems to be overlooked in the previous investigations in this field, 
in the fully non-perturbative treatment. We propose a trick which enables one to incorporate 
the Gribov horizon directly into the self-consistent Schwinger-Dyson equation in the gauge-fixed 
Yang-Mills theory, using the Gribov-Zwanziger framework with the horizon term. We obtain the 
following results, irrespective of the choice of the horizon term, (i) We find that there exists a fam- 
ily of solutions parameterized by one-parameter wr (0) which was assumed to be zero implicitly. 
The family includes both the scaling and decoupling solutions, and specification of the parameter 
discriminates between them, (ii) We observe that the inclusion of the horizon term cancels the 
ultraviolet divergence in the Schwinger-Dyson equation for the ghost propagator and the resulting 
(non-perturbative) self-consistent solution becomes ultraviolet finite. In other words, the horizon 
condition interpolates between the infrared behavior and the ultraviolet one. This leads to the in- 
compatibility of the conventional multiplicative renormalization scheme with the Gribov horizon 
condition. 
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1. Introduction 

We consider the SU(N) Yang-Mills theory in D-dimensional Euclidean space. A naive def- 
inition of a quantum Yang-Mills theory in the Landau gauge = is given by the functional 
integral: 

Z YM := J[d£f]8(d£f)det(-dD[£f]))exp{-S Y M[^]}, (1-1) 

where Sym is the Yang -Mills action and —dD[gf] is the Faddeev-Popov operator. For this formula 
to be correct, each gauge orbit must intersect the gauge fixing hypersurface Y := {.s/;ds^ = 0} 
once and a unique representative from each gauge orbit can be chosen by imposing ds? = 0. 
However, it has been shown that there exist many representatives obtained as intersections of a 
gauge orbit with the gauge fixing hypersurface T, called the Gribov copies [{[]]. 

In order to avoid the Gribov copies, Gribov proposed to restrict the functional integral to 
the 1st Gribov region £2: 

Z Gri bov:= / [d^]8{d£f)det(-dD[£f]))exp{-S Y M[^}}, (1-2) 

JQ. 

where the Gribov region defined by 

a:={£/;d&? = 0, -dD[af] >0} CT (1.3) 

is a bounded and convex region including the origin {&/ = 0}. In fact, —d^D^ \stf = 0] = — d^d^ > 
0, i.e., {&/ = 0} £fl. The boundary of £2 is called the Gribov horizon: 

dQ,:={s?f;dstf = 0, -dD\af\ =0}. (1.4) 

Then he predicted that the resulting Green functions exhibit unexpected behavior in the deep in- 
frared (IR) region, which plays the essential role in confinement. The gluon 2-point function (full 
or complete propagator) 



Df v (k):=8 AB 



k^k v \ F(k 2 ) a k^k v 



k 2 k 2 k 2 k 2 



(a = 0), (1.5) 



and the ghost propagator 



G AB (£):=-S AB ^, (1.6) 



k 2 

with the free case: F(k 2 ) = 1 and G(k 2 ) = 1 behave in the deep infrared (IR) region k 2 < 1 as 



2 



where M is a constant with mass dimension called the Gribov mass. This power like behavior 
should be compared with the ultraviolet (UV) behavior with the logarithmic corrections given by 
one-loop resumed perturbation or one-loop renormalization group: 

t 2 V /. k 2 \ S 13 



F(f)~(in^ , ■ r=- Tr t = - M - a* 
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The Gribov prediction was investigated more elaborately by solving the coupled Schwinger- 
Dyson (SD) equation for the gluon and ghost propagators where the dressing functions F and G 
are assumed to obey the power behavior with exponents a and /3 respectively: 

F(k 2 )=Ax(k 2 ) a , G{k 2 ) = B x (k 2 f , a + 2/3=0, 0<A,B<oo. (1.9) 

It has been shown that the scaling relation a + 2/3 =0 holds and hence a single exponent K is 
enough for characterizing the IR behavior: 

a = 2fc>l, j3 = -K<0, 1 /2 < K- < 1 (Gribov k= 1). (1.10) 

The Gribov prediction corresponds to the limiting value K = 1 . This result, gluon suppression and 
ghost enhancement in the IR region, leads to the running (ghost-antighost-gluon) coupling constant 
with a non- trivial IR fixed point: 

g 2 (k):=g 2 F{k 2 )G 2 {k 2 )^{0<)g 2 AB 2 {<oo) (k 2 ^0). (1.11) 

This solution is called the scaling solution. The gluon propagator F {k 2 ) jk 2 vanishes in the IR limit 
k 2 I 0, while the ghost propagator becomes more singular than the free case in the IR region, or 
the ghost dressing function G(k 2 ) diverges, i.e., 

G(0)=°°. (1.12) 

See the excellent review by Alkofer and von Smekal ^ for details. 

This IR behavior was considered to be reasonable from the viewpoint of color confinement. 
Due to Kugo and Ojima [Qj, all color non-singlet objects can not be observed or confined, in other 
words, only color singlet objects are observed, if a criterion m(0) = — 1 is satisfied in the Lorentz 
covariant gauge (a sufficient condition for color confinement). 1 It is claimed in [Q] that in the 
Landau gauge, the Kugo-Ojima criterion for color confinement u(0) = —1 is equivalent to the 
divergent ghost dressing function G(0) = °° , since in the Landau gauge G(0) is related to w(0) 
according to 

G(0) = [l + «(0)]- 1 . (1.13) 

Until 2006, it seemed that the scaling solution has been confirmed by the SD equation, the func- 
tional renormalization group equation and numerical simulations on lattice. This lead to the ghost 
dominance picture for color confinement. This strategy is called the Gribov-Zwanziger/Kugo- 
Ojima color confinement scenario. 

So far so good. However, these results are questioned by the Orsay group at Universite de Paris 
Sud. By careful analyses of the SD equation, so-called the decoupling solution was discovered [g|: 

F(k 2 )=A' x(k 2 f, G(k 2 ) = B' x (k 2 f , a' = l, p' = 0, 0<A , ,B , <oo, (1.14) 

which leads to the running coupling going to zero in the IR limit: 

g 2 (k):=g 2 F(k 2 )G 2 (k 2 )~g 2 A>B' 2 k 2 ^0 (k 2 ^0), (1.15) 

'Note that the Kugo-Ojima theory for color confinement does not take into account the Gribov problem and is based 
on the usual BRST formulation where the exact color symmetry and the well-defined BRST charge are assumed. 
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although its possibility was mentioned also in [Q]. Moreover, reexaminations of numerical simula- 
tions on large lattices [0, g, 10, 11, |T^, [jj], functional renormalization group equation [|l4|] and 



other methods [ |17| , |18| ] seem to support this result. 77ie decoupling solution implies that the gluon 
propagator goes to the non-zero and finite constant in the IR limit, while the ghost propagator 
behaves like free, namely, the ghost dressing function G(0) is non-zero and finite in the IR limit: 



< G(0) < °o. 



(1.16) 



In the decoupling solution, the gluon decouples below its mass scale and the ghost is still dominant, 
although the ghost dominance in the decoupling solution is weaker than that in the scaling solution. 



2. Exact relationship between the ghost dressing function and the Kugo-Ojima 
parameter 



We wish to point out that the relation ( |1.13[ ) used in the above discussion is not precise, since 
the exact relationship between the ghost dressing function and the Kugo-Ojima parameter is [23, 
2§ 

G(0) = [l + u(0) + w(0)]- 1 . (2.1) 
Therefore, the Kugo-Ojima criterion h(0) = — 1 is equivalent to the scaling solution G(0) = °°, 



if and only if w(0) = 0. In fact, the relationship (2.1) is a special case of the relationship between 
the ghost dressing function G{k 2 ) and the Kugo-Ojima function u(k 2 ) with an additional function 

w(k 2 y. 

G l (k 2 ) = 1 + u(k 2 ) + w(k 2 ). (2.2) 



This identity was first derived in [J4p and confirmed in fl23j , |24| ], while it was shown also in [|2£ 
based on a different method. Here two functions u and w are defined from the modified 1 -particle 
irreducible (mlPI) part as 



B\m\PI 



8^ v u(k ) + -j2~w(k 



?AB 



where u{k 2 ) agrees with the Kugo-Ojima function usually defined by 2 

((D^) A (g^ v x V) B ) k := (V - ^\ 8 AB u{k 2 ). 
The mlPI part is defined from the two-point function of the composite operators (See Fig. |I|) 

\A/„^ v . c2\B\ lAB/,\ i \AB, 



where 



A AB (k) :=((g^ x ^) A ^ c )f\^ D ) k {^ D {g^ v x tf) B ) k pl . 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



"This definition could be reconsidered from the viewpoint of renormalizability. 
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Figure 1: Diagrammatic representation of ((g-afjj x ^) A {g^ v x ((g^ x ^) A {g^ v x <^) B )|= onn , 

x x and x ?)^(^ v x PI . 



We should consider the renormalized version which will be discussed later: 

G R \k 2 ) = \ + u R (k 2 )+w R (k 2 ). 



(2.7) 



Here we understand that the renormalized quantities Gr, ur, and wr are free from the ultravio- 
let divergence. Note that wr(0) ^ means the existence of the massless pole in the correlation 
function X^ v (k). Indeed, wr(0) = is true at the tree level and holds in perturbation theory, since 
the massless pole coming from the elementary Faddeev-Popov ghost was already removed in the 
definition of the 1PI function X^ v {k), see Fig. [TJ Therefore, within the perturbation theory (up to a 
finite order), the scaling solution can be compatible with the Kugo-Ojima criterion: 



G s (0)=oo, UR (0) = -1, (w*(0)=0). 



(2.8) 



However, the Kugo-Ojima criterion ur (0) = — 1 can not be realized in perturbation theory, since 
it is a quantity of order g . Therefore, we must reexamine the whole strategy of GZ/KO from a 
fully non-perturbative point of view. For wr (0) ^ to be realized, a massless boundstate must be 
formed anew in the channel by the non-perturbative effect and the massless pole must be generated. 
However, the existence of such a pole would not mean the appearance of the physical massless 
particle in the spectrum, since it is not a gauge-invariant object. Therefore, wr(0) / does not 
contradict with the experiment. This consideration leads to a possibility of the decoupling solution 
satisfying the Kugo-Ojima criterion for color confinement. 

0<G«(0)<oo, UR (0) = -1, (w*(0)^0). (2.9) 

A new possibility of wr(0) ^ should be examined from a non-perturbative point of view, as 



pointed out in [25]. We can consider other possible cases which must be compatible with (2.7). 



3. The horizon condition and horizon terms 

In order to see the effect of the existence of the Gribov horizon, we make use of the Gribov- 
Zwanziger theory [19, 20, [2~I| ]: 



Z GZ := J [^]S(d^)det(^[^])exp J -S Y mW\ -J J d D xh(x) 1 , (3.1) 



where K is the Faddeev-Popov operator K[,e/] := — d^D^s^] = —d^(d^ +g^x) and h{x) = 
h[s/](x) is the Zwanziger horizon function. Here the parameter y called the Gribov parameter is 
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determined by solving a gap equation, commonly called the horizon condition: for D-dimensional 
Euclidean SU(N) Yang-Mills theory, 

(h(x)) GZ = (N 2 -\)D. (3.2) 

The horizon function plays the role of restricting the integration region inside the Gribov horizon 
d£l, since the horizon function has the property: 



£/ -, da K\sf\ | J d D xh[£/](x) T oo =^ exp{-7 J d D xh{x)} [ (y > 0). (3.3) 

However, the explicit form of the horizon function respecting exactly the original definition is not 
known and the choice of the horizon function is not unique at present. The first proposal of the 
horizon function is fli9| ] 

h(x) = J d D ygf ABC ^(x)(K- l ) CE (x,y)gf AFE ^ F (y). (3.4) 



The second proposal is Q21| ] 



h(x) = J d D yD^] AC (x)(K- l ) CE {x,y)D^] AE {y). (3.5) 

In any case, inclusion of the horizon term makes the theory non-local. 

It should be remarked that the horizon term does not uniquely fix the gauge, since there are 
still Gribov copies inside the 1st Gribov region. Recently, a one-parameter family of correlation 
functions are constructed in lattice gauge theory distinguished by a second gauge parameter B 



(Landau-B gauge) [16]. This uniquely specifies a representative from a gauge orbit and no further 



freedom in choosing a Gribov copy. 

4. Schwinger-Dyson equation with the horizon condition 

The SD equations for Green functions do not change their form even in the presence of the Gri- 
bov horizon, since the integrand of the functional integration formula for Green functions vanishes 
at the Gribov horizon as the boundary of the functional integration region due to the Faddeev- 
Popov operator. Hence, the solutions of the SD equation include both the solution with the Gribov 
horizon and the solution without restriction. Remarkably, it has been shown [|14|, [l5|] that the set 
of solutions of the coupled SD equation for gluon and ghost propagators is uniquely determined 
once a boundary value G(0) is given, corresponding to the scaling solution for G(0) = °° and the 
decoupling solution for < G(0) < °°. However, it is not yet examined how these solutions are 
related to the Gribov horizon. 

The Schwinger-Dyson (SD) equation for the ghost propagator i^€ Ac € B )^ in momentum space 



is written in the following form. We follow the notation of [|23j, g4j . 



which is obtained as the Fourier transform of 



<<rV~\- 1 = -8 AB k 2 - i^((g^ x tf) A <i B )f\ (4.i) 



= -((d^D^]^) A (x)^ B (y)} + 8 AB 8 D (x- y). (4.2) 
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In the Gribov theory, this is derived from the identity: 

8 



0= / \ d ^\[dm\w\w\-^y^\y) 

Jo. b c e A {x) L 



[ 8tf A {x) 

while in the Gribov-Zwanziger theory, the same form is obtained from 

0= / [djrt[dMdtf\[d&\ - 5 ir - \e-^-^ dDxh{x ^ B {y) 
J o^ A (x) L 



(4.3) 



(4.4) 



where 



Sym := ^ym + ^GF+FP, 



S Y u:=Jd D x^ v -^ v , 5 GF+FP := Jd D x{^-d^ + i^-d^D^}, (4.5) 
and the dot and the cross are defined as st ■ 3§ := $? A g3 A and (st x SS) A := f ABC $i B 3B c . 



In [25], we have given a trick to incorporate the horizon condition into the SD equation of the 
ghost propagator, which enables us to distinguish the solution associated with the particular choice 
of the horizon term and to discriminate the scaling and decoupling solutions. Following the idea 
of Gribov [[I]], we incorporate the horizon condition (h(0)) = (N 2 — l)D into the SD equation. By 
substituting the horizon condition into the free (tree) part of the SD equation, the SD equation for 
the ghost dressing function reads 

G- l (lc 2 ) = j^§ j5 + u(k 2 ) + W (k% (4.6) 

which is equivalent to the SD equation for the ghost propagator: 

<^V"\- 1 = -^J^y - §((8^ x tfjV)™. (4.7) 
Advantages of the SD equation with the insertion of the horizon condition is as follows. 

1. In order to obtain the scaling solution, the constant terms must cancel exactly or disappear 
at k = on the right-hand side of the SD equation. This is what implicitly assumed, but not 
stated explicitly, as pointed out by [^j]. Otherwise, the decoupling solution is obtained. The 
SD equation (^) with a horizon condition can be used to discriminate between the scaling 
and the decoupling solutions as a result of the existence of the Gribov horizon. 

2. The renormalized version of the SD equation still contains the UV divergence coming from 
the integration over the internal momenta in the self-energy part even after Green functions 
(the ghost propagator, gluon propagator and the vertex functions) are replaced by the renor- 
malized ones. If we include the horizon term, this UV divergence cancels between the self- 
energy term and the horizon term and the SD equation becomes a self-consistent equation to 
yield finite Green functions. 

We show that both horizon terms allow the existence of one-parameter family of solutions 
parameterized by a real number wr(0) which has been assumed implicitly to be zero wr(0) = 
in the previous investigations. Therefore, wr(0) plays the role of an additional non-perturbative 
gauge parameter which uniquely specifies the solution. We consider both the unrenormalized and 
renormalized versions of the proposed SD equation with the horizon condition being included. 
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5. Ultraviolet finiteness of the self-consistent solution of the Schwinger-Dyson 
equation with the horizon condition 

5.1 The first horizon term 



The first horizon term (|3.4|) yields [ 23 , 24] 



(h(0)) = - lim((^ x ^) A (g^ x tf) A ) k 

A:— >0 

= -{N 2 - 1) {Du(0) + w(0) - G(0)[«(0) + w(0)] 2 } , (5.1) 

and 

«°» - „ (0) _^(0) + G-'(0)-2 + G(0) (52) 



(N 2 -l)D D D 

Then the SD equation reads 

G - 1( , 2) _ O-'(0)-2 + C(0) _ ^ _ „ (0) + „ ( , 2) + „, (t2) . (5 3) 

Then, it is observed that the term — w(0) coming from the horizon condition exactly cancels the 
self-energy term u(k 2 ) at k = in the right-hand side of the SD equation. In the deep IR limit 
k = 0, therefore, we have 

G-'(0) = G ~' (0) - 2+g( V 4 + lW (5-4) 



D V D 

By solving this equation for G(0): G 2 (0) — [2 — (D — l)w(0)]G(0) + 1 — D = 0, we have 



G(0) = 1 - (D - l)w(0) /2 + V [1 - (D - l)w(0) /2] 2 - 1 + D > 0, (5.5) 
andM(O) = -l-w(0) + G- 1 (0) reads 



M (0) = -1 - w(0) - i 1 2 - 3w(0) - yjl2 + [2 - 3w(0)] 2 | . (5.6) 

This implies that the horizon condition determines the boundary value G(0) in the ghost SD equa- 
tion. Consequently, we have one-parameter family of solutions parameterized by w(0). 

We consider G(0) and u(0) as functions of w(0). Both G(0) and m(0) are monotonically 
decreasing functions in w(0); G(0),m(0) — > +o° as w(0) —> while G(0) — > and «(0) — » —5/3 
as w(0) — > +oo. Thus the scaling solution G(0) = +°° is obtained only when w(0) = — °°. Otherwise 
w(0) > the decoupling solution < G(0) < °° is obtained. 

Using a special value as an additional input w(0) = which is assumed in 23] or suggested 



by an independent analysis [E9], G(0) is determined self consistently by solving the above SD 



equation as [23, 24] 



G(0) = 1 + VD >0, u(0) = (-D±Vd)/{D-\). (5.7) 
In particular, for D = 4, 

G(0)=3>0, u(0)=-2/3 (D = 4). (5.8) 
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This agrees with the current data of numerical simulations [[31], 32]. 

We proceed to the treatment of the ultraviolet divergence. The naive SD equation without the 
horizon condition reads 

G A \k 2 ) = \ + u A (k 2 ) + w A (k 2 ). (5.9) 

Here, in order to avoid the ultraviolet divergence in the self-energy part u(k 2 ) + w(k 2 ), the UV 
cutoff A has been introduced thereby to make the self-energy part u A {k 2 ) -\-w\{k 2 ) finite. Conse- 
quently, G must depend on A, i.e., G A (k 2 ) := G(k 2 ,A). In this sense, the above SD equation is an 
unrenormalized version. 

It is pointed out in that by writing w(0) and w(0) as function of (7(0), one finds m a (0) — > 
+oo and wa(0) — > — °° such that ma(0) +wa(0) — > — 1, provided that G A (0) — > °° as A — > °°, see 
Appendix C of [^5]]. This analysis uses the first horizon condition and the relation G(0)~ l = 1 + 
m(0) + w(0). However, the statement G A (0) —> °° as A — ► oo is a result of perturbation theory. This 
analysis does not use the full information coming from the relation G(k 2 )~ l = 1 +u(k 2 ) +w(k 2 ) 
(the Schwinger-Dyson equation for the ghost dressing function) for the whole momentum region. 

In what follows, we show that G A (0) remains finite in a non-perturbative way, even after 
sending the UV cutoff to infinity A — » oo, once the full information on the whole momentum 
region is used. We consider the SD equation with the horizon condition: 

G A \k 2 ) = <r A \0)-2 + G A (0) _ w^O) _ ma(0) + uA{k2) +wA{](2y (5 1Q) 

If the horizon condition is incorporated into the SD equation, a partial cancellation at k = 
occurs between the horizon condition and the ghost self-energy. This cancellation always occurs 
for the (multiplicative renormalizable) part coming from A^)(0), while this is not the case for the 
contribution from the remaining term A^(0). 

The first horizon condition which is non-linear in G 1 is not multiplicatively renormalizable 



Q25Q. Therefore, the ultraviolet divergence can not be removed within the usual renormalization 
scheme. However, a novel situation occurs by introducing the first horizon condition into the SD 
equation. By the resulting combination u A {k 2 ) — w A (0), the ultraviolet divergence cancels exactly 
and the ultraviolet cutoff A can be sent to infinity to obtain a finite function of k 2 . Then the 
above SD equation is regarded as a self-consistent equation to give a finite ghost function G{k 2 ) = 
limA^°oGA(& 2 ) < °°- The term w\{k ) is finite from the beginning by some reason, although the 



A dependence is apparently assumed, see Appendix B of [25]. In other words, the SD equation 
is self-organized (in a non-perturbative way) to give a finite result. Therefore there is no need for 
the specific ultraviolet renormalization in the presence of the first horizon term! Thus, the results 
obtained in the unrenormalized case hold also after the ultraviolet cutoff A is send to infinity, as far 
as the ghost propagator or the dressing function is concerned. 

Even if the ghost dressing function is ultraviolet finite, we have still finite renormalization 
coming from the choice of the renormalization point, leading to consistent results with numerical 



simulations. See [25]. 



5.2 The second horizon term 



We can perform the similar analysis for the second horizon term. See [25] for details. 
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6. Conclusion and discussion 

We have discussed how the existence of the Gribov horizon affects the deep infrared behavior 
of the ghost propagator in the Landau gauge G = SU (N) Yang-Mills theory, using the Gribov- 
Zwanziger framework with the horizon condition (h(x)) = (dimG)D. Moreover, we have shown 
how to incorporate the horizon condition into the Schwinger-Dyson equation for the ghost propa- 
gator to discriminate between scaling and decoupling. We have examined two horizon conditions 



derived from two types of horizon terms, both of which were proposed by Zwanziger [ |19| , [21| ]. 
We have shown that one parameter family of solutions parameterized by w(0) exists in both cases, 
although some results crucially depend on the choice of the horizon term adopted. The value w(0) 
has been assumed implicitly to be zero w(0) =0 in previous studies. 

Thus the investigation of wr(0) is crucial to see which solution is realized. An interesting step 
towards this direction was done in p9[|. However, according to our analysis, their result, namely, 
the decoupling solution with wr(0) = is not compatible with the multiplicative renormalization 
scheme (without the horizon term): 

G-\k 2 )=Z c l G R \k 2 ) ^l + u(k 2 )=Z c l [l + u R (k 2 )], w(k 2 )=Z c l w R (k 2 ), (6.1) 

since w(k 2 ) is ultraviolet finite in contradict to Zq = °°. This issue should be reexamined and 
confirmed by further investigations. 

We can regard that wr(0) plays the role of an additional non-perturbative gauge parameter 
which uniquely specifies the solution from a one-parameter family of solutions including the scal- 
ing and decoupling. In view of this, the value of wr(0) itself might have no physical meaning. In 



fact, it has been shown [ |37| ] that all solutions (decoupling as well as scaling) lead to quark confine- 
ment by proving the vanishing of the Polyakov loop as a gauge-invariant order parameter of quark 
confinement. In this sense, discriminating between decoupling and scaling may not be so important 
from the physical point of view and main contribution to phenomenological studies comes from the 
high-momentum region and the intermediate region around lGeV which is stable irrespective of 
adopting the scaling or decoupling solution. It is shown that both scaling and decoupling solutions 



do not contradict the general principles of quantum gauge field theories [38]. 

How the existence of the horizon is relevant for color confinement. In the Gribov-Zwanziger 
theory (restricted to the 1st Gribov region), the BRST symmetry is broken by the existence of the 
horizon. 8Sgz = SSy ^ 0. Nevertheless, there exists a "BRST" like symmetry (without nilpotency 



39|] or with nilpotency [ |40| ] which leaves the Gribov-Zwanziger action invariant. Then we could 
apply the Kugo-Ojima idea to the Gribov-Zwanziger theory, which opens the path to searching 
for the modified color confinement criterion a la Kugo and Ojima. On the other hand, defining a 
non-perturbative BRST transformation without restricting to the first Gribov region will be another 



interesting possibility to be investigated [41] 
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